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Abstract. We consider four two- level atoms interacting with independent non- 
Markovian reservoirs with detuning. We mainly investigate the effects of the detuning 
and the length of the reservoir correlation time on the decoherence dynamics of the 
multipartite entanglement. We find that the time evolution of the entanglement of 
atomic and reservoir subsystems is determined by a parameter, which is a function 
of the detuning and the reservoir correlation time. We also find that the decay and 
revival of the entanglement of the atomic and reservoir subsystems are closely related 
to the sign of the decay rate. We also show that the cluster state is the most robust to 
decoherence comparing with Dicke, GHZ, and W states for this decoherence channel. 
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1. Introduction 

Entanglement plays a central role in quantum information processing. Since the 
maintenance of entangled states is very important in quantum information processing 
systems, the study of the effect of decoherence on entangled states would be of 
considerable importance from theoretical as well as experimental point of view. Recently 
entanglement of open quantum systems has attracted considerable attention due to its 
significance for both fundamentals and applications of quantum information processing 
[H |2l [3l m ini [HI [7] . Yu and Eberly [3] found that the decay of a single qubit coherence can 
be slower than the decay of entanglement, and they named the abrupt disappearance 
of entanglement at a finite time entanglement sudden death (ESD), which has been 
experimentally tested very recently [8]. Later Lopez et al. [5] revealed that when 
the bipartite entanglement suddenly disappears, the entanglement of corresponding 
reservoir suddenly and necessarily appears, which is called entanglement sudden birth 
(ESB). Although the bipartite entanglement is well understood in many aspects, 
until now the multipartite entanglement is far from clear and thus deserves profound 
understandings. Several important applications of multipartite entangled states, like 
quantum error correction, quantum computing, etc. have been found (for recent reviews 
see Refs.[9l[T0l|ll] and references therein). Recently, people investigated the decoherence 
dynamics of multipartite entangled states, for instance, Giihne, Bodoky and Blaauboer 
[12] discussed some multipartite entanglement properties under the influence of a global 
dephasing process using the geometric measure of entanglement, and they showed that 
the Dicke state is the most robust to decoherence comparing with GHZ, cluster, and W 
states. Borras and Majtey et al. [13J investigated the robustness of highly entangled 
multiqubit states under different decoherence channels [llj, and later, they further 
studied the geometry of robust entangled multiparticle states under decoherence [T5] . 

In realistic physical systems the assumption of Markovian dynamics can only 
be an approximation. Generally speaking the quantum system of interest interacts 
with the reservoirs, whose spectral density strongly varies with frequency, which 
is called non-Markovian open quantum systems [16]. Non-Markovian dynamics is 
characterized by the existence of a memory time scale during which information that 
has been transferred from the system to the environment can flow back into the 
system. The non-Markovian systems appear in many branches of physics, such as 
quantum optics [Ml [ITJ [18], solid state physics [19], and quantum chemistry [20] . 
In quantum information processing, the non-Markovian character of decoherence has 
to be considered [21]. The non-Markovian dynamics has drawn great attentions 
including continuous- variable [22l [231 EH [251 EH [271 [2H] and discrete-variable systems 
[29l[30l[3ll[32l[33l[32l35l[36l[37]. Very recently N qubits initially in the mixed GHZ- 
type and W-type states interacting with independent structured reservoirs have been 
investigated and it is found that the A^-qubit entanglement revivals are related to the 
qubit number and the purity of the initial state of A^ qubits [35]. The measure of the 
degree of non-Markovian behavior [38] and a basic relation between the quantum Fisher 
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information flow and non-Markovianity have been proposed for quantum dynamics of 
open systems [39] . 

In this paper, we are very interested in the decoherence dynamics of multipartite 
entanglement under the influence of non-Markovian reservoir. We mainly investigate the 
effects of the detuning and the length of the reservoir correlation time on the decoherence 
dynamics of the multipartite entanglement. We consider that four noninteracting two- 
level atoms interact with four independent non-Markovian reservoirs with detuning. We 
let the four atoms initially be prepared in different four particle entangled states (cluster, 
Dicke, GHZ, and W states). We find that the dynamical behaviors of the entanglement 
of the atomic and reservoir subsystems for various initial states are determined by a 
complex parameter, which is a function of the detuning and the reservoir correlation 
time. The real part of this parameter is closely related to the decay of the entanglement 
and the imaginary part of it is closely related the oscillations of the revivals, and all the 
dynamical behavior of the entanglement can be uniformly explained by this parameter. 
We also find a close relationship between the decay rate and the entanglement of the 
atomic and reservoir subsystems, i.e., the decay or revival of the entanglement of the 
atomic and reservoir subsystems are decided by the sign of the decay rate. We also 
find that the cluster state is the most robust to decoherence, which is different from 
the results of Ref . [T2] . Here we consider the amplitude damping process, while Ref.[T2] 
discuss the global dephasing process. The paper is organized as follows: In Sec. 2, 
we present the model which consists of four atoms interacting with individual non- 
Markovian reservoirs, and introduce the multipartite entanglement measure. We present 
and analyze the results in Sec. 3. Finally, we give some conclusions in Sec. 4. 

2. Model 

We consider a system of four identical noninteracting two-level atoms, each of them 
coupled to its own reservoir. For simplicity, we assume that each corresponding reservoir 
is initially in the vacuum state. Due to the independence of each atom, we only need to 
discuss the problem of a single atom interacting with its corresponding reservoir. The 
Hamiltonian of the interaction between a single atom and A'^-mode reservoir under the 
rotating- wave approximation can be written as {h = 1), 

N N 

H = UaCT+a^ + ^ ujkaldk + ^ 5'fc('5"-4 + dr+Ok), (1) 

k=l k=l 

where (T+ = |1)(0| and o"_ = |0)(1|, are the Pauli raising and lowering operators for the 
atom respectively, and Ua is the transition frequency of each atom. a|. and dk are the 
creation and annihilation operators with frequency Uk for the reservoir mode k, and gk 
is the corresponding coupling constant. 

Here we consider one excitation case, namely, the atom and corresponding reservoir 
are initially in the excited state and vacuum state respectively, i.e., \ipo) = \l)a ® \0)r, 
where |0),. = H^i |OA;)r- The subscripts a and r refer to the atom and the corresponding 
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reservoir respectively. So the state of the total system at any time t can be denoted by 

N 

\iJt) = U{t)\l)a\0)r + Y,Dk{t)\0)a\lk)r, (2) 
fc=l 

where the state \lk)r represents the reservoir having one excitation in mode k. 

Similar to the method used in Ref. [16], we can obtain a closed equation for the 
coefficient z/(t) in Eq.(2), 

t 

m = - j f{t-h)y{ti)dt^, (3) 



where the kernel fit — ti) = f duJ{ijj)exp[i{ijJa — ou)(t — ti)] is related to the spectral 
density of the reservoir. We take a Lorentzian spectral density of the reservoir 

mm 



the center of which is detuned from the transition frequency Ua of the two-level atom 
by an amount 6. And the parameter A defines the spectral width of the coupling, 
which is associated with the reservoir correlation time by the relation tb = and the 
parameter 70 is related to the relaxation time scale tr by the relation tr = Jq^. 

By making the Laplace transformation of Eq.(3), we can obtain the solution of 
with initial condition z/(0) = 1, 

u{t) = e-(^-*'^)*/2 [cosh(xt/2) + —— sinh(xt/2)] , (5) 

where 

X = V(A-^<5)2-27oA. (6) 

Furthermore, if we take the form of a collective state of the reservoir, namely, letting 
= {l/^i(t)) J2k=i Dk{t)\l)r with = - Eq.(2) can be rewritten as 

\ijt) = l^{t)\l)a\0)r + Ktma\l)r. (7) 

There are many entanglement measures to quantify the bipartite entanglement such 
as von Neumann entropy, negativity [lOl HI] , concurrence [12] , relative entropy [13] , etc. 
However, generally speaking, for multipartite system the definition of entanglement 
measure is difficult. Up to now, there are some entanglement measures proposed 
to quantify the multipartite entanglement, which include Schmidt measure [S], the 
geometric measure of entanglement [15], the global entanglement measure [l6l W7\. etc. 
In this paper, we will use a very popular entanglement measure [131 115] which is averaged 
over all possible bipartitions. The mathematical definition of this measure is 

^= ^-L^ ^E^'"), (8) 



[N/2] 



m=l 
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where 

^ bipart 

bipart 

E{i) indicates the entanglement connected to single bipartition of the A^-qubit system, 
and E^'^''^ denotes the average entanglement over all nonequivalent bipartitions iVbipart 
between subsets of m qubits and the remaining N — m qubits. If one uses the linear 
entropy Sl of the reduced density matrix of the smaller bipartitions to compute E{i), it 
will reduce to the well known Meyer- Wallach multipartite entanglement measure 
Because in this paper we will deal with mixed states, we take the negativity to measure 
the bipartite entanglement. The normalized negativity is defined as |T5] 



i 

where is the negative eigenvalue of the partial transpose matrix for m and the 
remaining N — m bipartition. 

3. Results and discussions 

We suppose that four two-level atoms are initially prepared in some famous multipartite 
entangled states (cluster, Dicke, GHZ, and W states), and the corresponding reservoirs 
are initially prepared in the vacuum state. For W state the initial state of atom-reservoir 
system is: 

|0o) = (|0001)„,,,,3,, + |0010),,„,„3„, + |0100)„,,,,3,, 

+ |1000)„,„,„3„J|0000),,,,,3,,/2, (11) 

where the subscripts ai{i = 1,2,3,4) represents the atom, and ri{i = 1,2,3,4) refers 
to the corresponding reservoir. From Eq.(7), the evolution of the total system can be 
obtained 



10*) = [|000),^„^„3 1000),^,^,^ (z/(t)|l),jO)., +/i(t)|0),Jl), 

+ Wa,a,a, W r.r^r, (K^) I l)aj0),J + /i(t) |0)„, 1 1),J 

+ |000),^,3,^ |000),^,3,^ (i/(t)|l),jO),, +/i(t)|0),Jl),J]/2. (12) 

So the reduced density operator of atomic subsystem Pait) = Tr^d^i) (0t|) and the 
reduced density operator of reservoir subsystem pr(t) = Tr a{\(pt) {4>t\) can be obtained. 
Then from Eq.(8-10) we can obtain the degree of the entanglement of the atomic and the 
reservoir subsystems. For convenience, we denote the degree of the entanglement of the 
atomic subsystem and the corresponding reservoir subsystem by Ea and Er, respectively. 
Here for W state, Ea and Er can be obtained, 

Ea = |[16-16|z/(t)|'-(6V7|z/(t)|^ - 8\u{t)\^ + 4+4^21^^)1^ - 2\u{t)\^ + l)]|/24, (13) 
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= |[16|z/(t)|'-(6V7|Kt)|' - 6|z/(i)P + M^2\u{t)\^ - 2|z/(t)|2 + 1)]|/24.(14) 

Similarly we can obtain Ea and Er for the other initial states, cluster state ICL4) = 
(|0000) + |0011) + |1100) - |llll))/2 , Dicke state ID4) = (|0011) + |0101) + |1001) + 
|1100) + |0110) + |1010))/V6 , and GHZ state \DHZi) = (|0000) + |1111))/V2, 
respectively. 

First we consider the resonant case, i.e., 5 = 0. (i) A = IO70, which corresponds 
to the Markovian regime. In Fig.l we plot the entanglement evolution of the atomic 
and reservoir subsystems for initial cluster, Dicke, GHZ, and W states. From Fig. la 
we can see that the cluster state is the most robust against decoherence. The result is 
quite different from that of Ref . [12] , in which the Dicke state is the most robust against 
decoherence among GHZ, cluster, and W states. This is because that the decoherence 
channel used in Ref. [12] is different from ours. In Ref. [12] the global dephasing channel 
is considered, while in this paper we consider the amplitude damping channel. This 
means that the most robust multipartite entangled state might be different for different 
decoherence channels. It should be noted that the entanglement measure different 
from ours is used in Ref. [T2] . but we use our entanglement measure to recalculate the 
entanglement of Ref. [12] , and the result is the same as that of Ref. [12] . From Fig.l we can 
see that the entanglement of the atomic subsystem for all the initial states decreases 
monotonically to zero, the entanglement of the reservoir subsystem for all the initial 
states increases monotonically to the steady maximum, and the entanglement contained 
initially in the atomic subsystem is finally transferred into the reservoir subsystem. 

(ii) A = O.I70, which corresponds to the non-Markovian regime with relatively 
short reservoir correlation time. In Fig. 2 we plot the entanglement evolution of the 
atomic and reservoir subsystems for initial cluster, Dicke, GHZ, and W states with 
A = O.I70 and 6 = 0. The dynamical behaviors of the entanglement of the atomic and 
reservoir subsystems in the non-Markovian regime are quite different from that in the 
Markovian regime. It can be seen from Fig. la and Fig. 2a that the common feature 
for different initial states in the non-Markovian regime is that the entanglement of the 
atomic subsystem decreases to zero much slowly than that in the Markovian regime, and 
in the non-Markovian regime after the entanglement of the atomic subsystem decays to 
zero it can revive at later time, which is quite different from the Markovian case. The 
reason is that the information, which the atomic subsystem loses to the reservoir, is 
later recovered by the atomic subsystem due to the reservoir non-Markovian memory. 
It is noted that in both the Markovian and non-Markovian regimes all the initial 
entanglement of atomic subsystem Ea will decay and is eventually lost for long times, and 
the entanglement of reservoir subsystem Er gradually increases to the steady maximum 
from zero, which can be seen from Figs.l and 2. It can be seen from Figs, lb and 2b that 
in the non-Markovian regime the entanglement of the reservoir subsystem Er at first 
shows oscillations as a function of time for all the initial atomic states, and finally the 
steady maximum entanglement is achieved, while in the Markovian regime Er increase 
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to the steady maximum monotonically. All the distinction between the entanglement 
properties in the Markovian regime and that in the non-Markovian regime is induced by 
the non-Markovian memory. In other words, in the Markovian regime the information 
flow is one directional, namely from atoms to reservoirs, while in the non-Markovian 
regime the information flow is bidirectional, namely the exchange of information back 
and forth between the atomic and reservoir subsystems, which causes the oscillations of 
the entanglement of the atomic and reservoir subsystems. 

(iii) A = O.OI70, which corresponds to the non-Markovian regime with the relatively 
long reservoir correlation time. In Fig. 3 we plot the entanglement evolutions of atomic 
and reservoir subsystems for the four initial states with A = O.OI70 and 5 = 0. 
Comparing Figs. 2a and 3a, we can find that the revival of the with relatively 
long reservoir correlation time is more obvious than that with relatively short reservoir 
correlation time, i.e., the amplitude of revival with relatively long reservoir correlation 
time is much larger than that with relatively short reservoir correlation time. For 
reservoir subsystem, compared with the case with relatively short reservoir correlation 
time, it is more difficult to achieve the steady maximum of entanglement with relatively 
long reservoir correlation time, which can be seen from Figs. 2b and 3b. This can be 
understood as follows: Increasing the reservoir correlation time means that the memory 
effect of the reservoir becomes stronger, and then the amount of information exchanged 
between the atomic and the reservoir subsystems will be enhanced. So the atomic 
subsystem can obtain more information from the reservoir subsystem in the case of 
relatively long reservoir correlation time and the revival is stronger, and because of the 
enhanced information exchange back and forth the reservoir subsystem will need more 
time to achieve the final maximum entanglement. 

Now we consider the off-resonant case, i.e., 5 — 8A. (i) A = O.I70. In Fig.4 we 
plot the entanglement evolution of atomic and reservoir subsystems for the four initial 
states with 5 = 8A and A = O.I70. In the off- resonant case the entanglement of the 
atomic subsystem Ea decays to zero with small oscillations, and during each oscillation 
Ea can not collapse to zero. And the overall decay rate is smaller than that in the 
corresponding resonant case, which can be seen from Figs. 4a and 2a. The entanglement 
of the reservoir subsystem E,. at first increases with very small amplitude oscillations in a 
very short period of time and then increases monotonically to the steady entanglement, 
and the overall increasing rate is smaller than that in the corresponding resonant case, 
which can be seen from Figs.4b and 2b. This can be easily understood: When the value 
of the detuning 5 increases, the effective coupling between the atomic and reservoir 
subsystem decreases. So the exchange of information between the atomic subsystem 
and the reservoir subsystem is not effective and adequate, (ii) A = O.OI70. In Fig. 5 we 
plot the entanglement evolution of atomic and reservoir subsystems for the four initial 
states with 5 = 8A and A = O.OI70. Comparing Figs. 5b and 3b we can find that due 
to the increasing of 5, the exchange of information is not effective, the oscillations of 
Er are not adequate, more specifically Er can not achieve its maximum during each 
oscillation. From Figs. 5a and 3a it can be found that increasing the detuning 5 the 
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period of the revival is shorten, and the amphtude of the revival increases. This result 
is very interesting. As we have mentioned above increasing the detuning 6 will make 
the exchange of information less effective, then why the revival becomes stronger? Now 
we analyze the decoherence dynamics of the multipartite entangled states in detail. 

To gain insight in the physical processes characterizing the decoherence dynamics 
for different initial states, we consider the parameter x, and find that all the above 
phenomenon can be uniformly explained by this parameter. From Figs. 1-5 it is easy 
to find that the dynamical behaviors of the entanglement for different initial states are 
very similar. For simplicity in the following we will take W state as an example. From 
Eq.(6) we can see that generally x is a complex number, and we will show that the real 
part Rex is responsible for the decay of Ea and the imaginary part Imx is responsible 
for the oscillations associated with the revival. From Eq.(13) we can see that the degree 
of entanglement is a function of which means that all the decoherence dynamics 

of the entanglement entirely depends on And it is easy to see from Eq.(5) that in 
the long time limit u{t) is dominated by the terms containing the factor e^"^"*''^'^^!^*/^. 
From numerical calculations we find that |Rex| increases with 6 and is always less than 
A, and |Imx| also increases with 6. From Eqs.(5) and (13) roughly speaking, A — |Rex| 
determines the decay of the entanglement, which we call it the decay exponent to be 
distinguished from the decay rate 7(t) [16], and |Imx| determines the basic frequency 
of the oscillations in the revivals (it is noted that the overall phase factor e*''*/^ in u(t) 
does not make any contributions to the entanglement). When A > 270 and 6 = 0, from 
Eq.(6) X is a real number, i.e., |Imx| = 0, which corresponds to the Markovian regime. 
Hence Ea will decay exponentially to zero without oscillations, and the revival can not 
appear. It is easy to prove that when A > 270, the decay exponent A — y^A^ — 270A is a 
decreasing function of A, and approaches 70 with the increasing A, the maximum value 
of which is 270 occurring at A = 270. When 6 = and A < 270, which is corresponding 
to the non-Markovian regime, x is a pure imaginary, and the oscillations appear. In this 
case the decay exponent is just A. That is why the entanglement decay for A = O.I70 
corresponding to the non-Markovian regime (Fig. 2a) is slower than that A = IO70 
corresponding to the Markovian regime (Fig. la). And also that is why the entanglement 
with A = O.OI70 and 6 = (Fig. 3a) decays much slowly than that with A = O.I70 and 
6 = (Fig. 2a). Remember that |Rex| and |Imx| increase with 6, so the decay exponent 
A — I Rex I decreases with the increasing of 6. In this way the envelope of Ea{t) decay 
more and more slowly with the increasing of 6, so that during each revival the amplitude 
achieved is increasing with the increasing of 6. This explains why with the increasing 
the detuning 6 the period of the revival is shorten, and the amplitude of the revival 
increases (see Figs. 3a and 5a). Now we consider the dispersive regime, i.e., 6 ^ A, 70, 
and in this case in the long time limit i/(t) ^ 1 — iX'^/46^ and the steady entanglement 
of the corresponding atomic subsystem Ea ~ [6-\/3 + 4 + (20 + 6a/3)A^/165^]/24 can 
be achieved. This means that in the dispersive regime the decay of entanglement Ea is 
strongly inhibited. We also calculate the degree of entanglement for the initial five (and 
six) particle W state, and we find that the entanglement dynamics of atomic subsystem 
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for the initial five (or six) particle W state is almost the same as that of the initial four 
particle W state, more specifically the influence of the detuning and the length of the 
reservoir correlation time on the dynamical behavior of the entanglement for the initial 
five (or six) particle W state is almost the same as that for the initial four particle W 
state. Whenever the degree of entanglement for the four particle case increases, the 
degree of entanglement for the corresponding five (or six) particle case also increases, 
and whenever the degree of entanglement for the four particle case decreases, the degree 
of entanglement for the five (or six) particle case also decreases. And the time, at which 
the entanglement reaches the maximum (or the minimum), is almost the same for all 
the three cases. This can be understood. Because we find that the entanglement of the 
atomic subsystem for five (or six) particle case is also a function of |z/(t)p, which means 
that the dynamical behavior of the entanglement for the five (or six) particle is also 
decided by the real part and imaginary part of 

It is well known that in the most general form of a time-local master equation for the 
reduced density operator, the decoherence is induced by the Lindblad (jump) operator 
with a decay rate 7(t). If the decay rate 7(t) is always positive, this describes the so- 
called time- dependent Markovian process [Ml EHl HQ] , but if at least during a period of 
time the decay rate 7(t) is negative, the non-Markovian process emerges. Now we also 
take W state as an example to show the relation between the decoherence dynamics of 
the entanglement and the decay rate 7(t). For simplicity we let A = O.OI70 and 5 = 0, 
and in this case the decay rate 7(t) can be expressed as 7(t) = — 2Re{z/(t)/t/(t)} [16] . 
where p{t) is obtained by choosing 5=0 in Eq.(5). In Fig. 6 we plot Ea, Er and j{t) as 
functions of scaled time 7ot for A = O.OI70 and 6 = 0. From Fig. 6, it is obvious to see that 
whenever ^{t) (dotted line) takes negative values, Ea (solid line) begins to revive and 
increase monophonically, and the corresponding E^. begins to decrease monophonically; 
when 7(t) takes positive values, Ea will begin to decrease monophonically, and Er begins 
to increase monophonically. This can be easily understood: When 7(t) is positive, the 
information flow is from atomic subsystem to reservoir subsystem, which means that 
Ea will decay, and E,. will increase; When 7(t) is negative corresponding to the memory 
effect of the reservoir, the information flow is from the reservoir subsystem to atomic 
subsystem, so Ea will revive and Er will decay. 

4. Conclusions 

In this paper, we have considered four atoms with initial entanglement interact with 
independent non-Markovian reservoirs. We have analyzed the decoherence dynamics 
for various initial states in Markovian (A = IO70), weak non-Markovian (A = O.I70) and 
strong non-Markovian (A = O.OI70) regimes, with and without the detunings. We have 
found that the decoherence dynamics of the atomic and reservoir subsystems strongly 
depends on a parameter, which is decided by the detuning and the reservoir correlation 
time, and all the phenomenon can be explained by this parameter. The real part of 
this parameter determines the decay the entanglement and the imaginary part of it 
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determines the oscillations of the revival. We also have found that whenever ^(t) takes 
negative values, Ea will begins to revive, and the corresponding Er begins to decrease; 
when 7(t) takes positive values, Ea will begin to decay, and Er will begin to increase. 
We have also found that for this decoherence channel the cluster state is the most robust 
to decoherence comparing with Dicke, GHZ, and W states. 
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CAPTIONS 

Figure 1: In the Markovian regime (A = IO70, 5 — Ea and as functions of 
scaled time ^^t for various initial states: (a) the atomic subsystem; (b) the reservoir 
subsystem. 

Figure 2: In the non-Markovian regime with relatively short reservoir correlation 
time (A = O.I70, 5 = Q) E^ and as functions of scaled time 70^ for various initial 
states: (a) the atomic subsystem; (b) the reservoir subsystem. 

Figure 3: In the non-Markovian regime with relatively long reservoir correlation 
time (A = O.OI70, 5 = 0) Ea and E,. as functions of scaled time 7ot for various initial 
states: (a) the atomic subsystem; (b) the reservoir subsystem. 

Figure 4: In the non-Markovian regime with relatively short reservoir correlation 
time with detuning (A = O.I70, 5 — 8A) Ea and Ej. as functions of scaled time 70^ for 
various initial states: (a) the atomic subsystem; (b) the reservoir subsystem. 

Figure 5: In the non-Markovian regime with relatively long reservoir correlation 
time with detuning (A = O.OI70, 5 = 8A) Ea and E,. as functions of scaled time ^ygt for 
various initial states: (a) the atomic subsystem; (b) the reservoir subsystem. 

Figure 6: The atomic entanglement Ea, reservoir entanglement Ej. and the decay 
rate 7(t) for initial W state as functions of •jot (A = O.OI70 and 6 = 0): the atomic 
subsystem (solid line); the reservoir subsystem (dashed line); the decay rate (dotted 
line) . 
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